Maier-Saupe-type theory of ferroelectric nanoparticles in nematic liquid crystals 
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Several experiments have reported that ferroelectric nanoparticles have drastic effects on nematic 
liquid crystals — increasing the isotropic-nematic transition temperature by about 5 K, and greatly 
increasing the sensitivity to applied electric fields. In a recent paper [L. M. Lopatina and J. V. 
Selinger, Phys. Rev. Lett. 102, 197802 (2009)], we modeled these effects through a Landau theory, 
based on coupled orientational order parameters for the liquid crystal and the nanoparticles. This 
model has one important limitation: Like all Landau theories, it involves an expansion of the 
free energy in powers of the order parameters, and hence it overestimates the order parameters 
that occur in the low-temperature phase. For that reason, we now develop a new Maier-Saupe- 
type model, which explicitly shows the low-temperature saturation of the order parameters. This 
model reduces to the Landau theory in the limit of high temperature or weak coupling, but shows 
different behavior in the opposite limit. We compare these calculations with experimental results 
on ferroelectric nanoparticles in liquid crystals. 



I. INTRODUCTION 

One important goal of modern liquid-crystal research 
is to enhance the properties of liquid crystals through 
physical methods, without the need for new chemical 
synthesis. One way to achieve this goal is to put col- 
loidal particles into liquid crystals. If the particles have 
a length scale of microns, they induce elastic distortions 
of the liquid crystals, and these distortions mediate an 
effective interaction between the particles. The particles 
may then form a periodic array, leading to a composite 
material with potential applications in photonics [H-Q- 
By comparison, if the particles have a length scale of 
10-100 nm, they are too small to distort the liquid crys- 
tal. In that case, the system enters into another regime 
of behavior, in which the particles function as molec- 
ular additives to change the effective properties of the 
liquid-crystal host. One particularly interesting case oc- 
curs if the particles are ferroelectric. Experiments have 
shown that low concentrations of ferroelectric Sn2P2Sg 
or BaTiOs nanoparticles increase the orientational or- 
der parameter, increase the isotropic-nematic transition 
temperature, and decrease the switching voltage for the 
Frederiks transition 04HI- Thus, they provide a new op- 
portunity to enhance the properties of liquid crystals for 
technological applications. 

To make further progress with these materials, it is es- 
sential to develop a theory for the interaction between liq- 
uid crystals and ferroelectric nanoparticles. In previous 
theoretical research, Reshetnyak et al. have developed a 
theoretical approach based on electrostatics [IH, [H, [l5[ . 
In this theory, the key issue is how an ensemble of 
nanoparticles with aligned dipole moments can polarize 
the liquid-crystal molecules, hence increasing the inter- 
molecular interaction. This electrostatic effect enhances 
the isotropic-nematic transition temperature and reduces 



the Frederiks transition voltage. In related research, 
Pereira et al. have performed molecular dynamics sim- 
ulations of ferroelectric nanoparticles immersed in a ne- 
matic liquid crystal [16]. These simulations also assume 
that the nanoparticles are aligned, and they also find a 
substantial enhancement of liquid-crystal order. 

In a recent paper [17j . we proposed a different type 
of theory for the statistical mechanics of ferroelectric 
nanoparticles in liquid crystals. In that theory, we sup- 
pose that both the liquid crystals and the nanoparti- 
cles have distributions of orientations, as illustrated in 
Fig. 1. These distributions are characterized by two ori- 
entational order parameters, which interact with each 
other. Using a Landau theory, we showed that the cou- 
pling stabilizes the nematic phase. By estimating the 
strength of the coupling, we calculated the enhancement 
in the isotropic-nematic transition temperature. We also 
predicted that the nanoparticles would greatly increase 
the Kerr effect, the response of the isotropic phase to an 
applied electric field. 

Although the work of Ref. [13] demonstrates an im- 
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FIG. 1: (Color online) Schematic illustration of ferroelectric 
nanoparticles suspended in a liquid crystal. The electrostatic 
dipole moments of the nanoparticles have a distribution of 
orientations. 
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portant physical mechanism, we must acknowledge that 
it has one mathematical limitation: Like all Landau theo- 
ries, it involves an expansion of the free energy in powers 
of the order parameters. This expansion is valid when the 
order parameters are small, but it breaks down when they 
become large. In particular, the theory allows the order 
parameters to become larger than 1 , which is clearly im- 
possible. For ferroelectric nanoparticles in a liquid crys- 
tal, the nanoparticle order parameter is not necessarily 
small, even near the isotropic-nematic transition. 

The purpose of the current paper is to generalize the 
previous theory by eliminating the assumption that the 
order parameters are small. For this generalization, we 
now use a Maier-Saupe-type theory instead of a Lan- 
dau theory. We still consider the same physical concept 
of coupled orientational order parameters for the liquid 
crystals and the nanoparticles, and we still use the same 
energy of interaction between them. However, we now 
use a more general expression for the entropy, not a power 
series, which enforces the constraint that the order pa- 
rameters cannot become larger than 1. This change al- 
lows us to avoid the potential mathematical inconsistency 
of Landau theory. 

Like our previous calculation, the work presented 
here shows that doping liquid crystals with ferroelectric 
nanoparticles enhances the isotropic-nematic transition 
temperature. In the limit of weak coupling between the 
nanoparticles and the liquid crystal, the Maier-Saupe- 
type theory exactly reduces to the Landau theory. How- 
ever, in the case of strong coupling, the new theory pre- 
dicts a smaller but still substantial enhancement. Rough 
estimates suggest that the experimental system is in the 
limit of strong coupling, so it is important to use this 
modified theory. Furthermore, the work presented here 
also predicts the Kerr effect as a function of applied elec- 
tric field. In the limit of low electric field, the Maier- 
Saupe-type theory exactly reduces to the Landau theory. 
However, for larger field, the nanoparticle order saturates 
and the enhanced Kerr effect is cut off. 

The plan of this paper is as follows. In Sec. II 
we present the formalism of Maier-Saupe theory, with 
interacting orientational distributions for liquid-crystal 
molecules and nanoparticles. In Sec. Ill we apply this 
formalism to calculate the isotropic-nematic transition 
temperature, and determine the enhancement due to 
nanoparticles. In Sec. IV we use the same formalism to 
calculate the Kerr effect of induced orientational order 
under an applied electric field, and investigate how this 
effect depends on the magnitude of the field. Finally, in 
Sec. V we discuss the main conclusions of this study. 



tal equation of mean-field theory, 

F= (H) + k B T(lng), 



(1) 



where the first term is the energy, the second term is the 
entropic contribution to the free energy, and the averages 
are taken over the single-particle distribution function g. 
Hence, the first step is to define the distribution functions 
for liquid-crystal molecules and nanoparticles. 

Liquid-crystal molecules are rod-shaped objects, with 
each molecule characterized by the direction of its long 
axis m. In the nematic phase, these axes are preferen- 
tially oriented along the average director n. Because the 
individual molecules are equally likely to point along +n 
or n, the single-molecule distribution function can be 
written as 



exp(C/ LC P2(cos0)) 



J\ rf(cos 6>) cxp([/ LC P 2 (cos 6>)) ' 



(2) 



where is the angle between the molecular orientation 
m and the average director n, and P2 is the second Leg- 
endre polynomial. The parameter {7lc is a variational 
parameter, which acts as an effective field on the molec- 
ular orientation. It is related to the standard nematic 
order parameter Slc = (P2{cos8)) by 

j\ d(cos 9)P 2 (cos 6>) cxp([/ LC P 2 (cos 0)) 
ShC ~ J\ d(cos 0) cxp(£/ LC P 2 (cos 0)) ^ 

Note that Ulc ranges from to 00, while Slc ranges 
from to 1. 

We can now calculate the free energy of a pure liquid- 
crystal system. Maier-Saupe theory assumes that the in- 
teraction energy between neighboring molecules i and j 
is proportional to —(mi • m.,-) 2 . With the assumed distri- 
bution function, the average interaction energy becomes 



^energetic - ( H ) = ~^J N hcSl C , 



(4) 



where -/Vlc is the number of liquid-crystal molecules in 
the system, and J is an energetic parameter proportional 
to the interaction strength and the number of neighbors 
per molecule. Furthermore, the entropic contribution to 
the free energy can be written as 



entropic 



k B TN LG (^g hC ) 



k B TN LC 



LC 



(5) 



In 



j* ^ d(cos 6) exp([/ LC P 2 (cos 0)) 



II. OVERVIEW OF MAIER-SAUPE THEORY 

In this section we introduce the free energy for a system 
of liquid-crystal molecules with ferroelectric nanoparti- 
cles. To construct the free energy, we use the fundamen- 



By combining these pieces, we obtain the total free en- 
ergy of the liquid crystal, 



J 



N LC k B T 



3k B T 



Slc + UlcSlc 



(6) 



In 



J\ d(cos0)exp([/ LC P 2 (cos0)) 
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The free energy of Eq. ([6]) is a function of the tem- 
perature T and the variational parameter L/lc > with Slc 
defined implicitly as a function of Ulc through Eq. ((3|). 
By minimizing the free energy over £/lc for varying tem- 
perature, we can find the liquid crystal has a first-order 
transition from the isotropic phase with Ulc = Slc = 
to the nematic phase with 

U LC = 1-95, S LC - 0.429. (7) 

The numerical solution for the transition temperature in 
this pure liquid crystal is 

J 



T N i = 0.147 



k-Q 



(8) 



Also, we can find an analytic solution for the limit of 
supercooling, 



2T T 
T* = J±- =0.133f- 
15kb 



(9) 



From experiments we know T* and Tni for any particular 
liquid-crystal material, so we can use Eq. © or © to 
determine J for that material, 



J = 6.81fc B Ti 



NI- 



(10) 



Once we add nanoparticles to the system, we get an- 
other distribution function for the orientations of the 
nanoparticle dipole moments. By symmetry, we expect 
that this distribution should be aligned along the same 
axis n as the liquid-crystal distribution. However, the 
magnitude of the order may be different. Hence, we can 
write the nanoparticle distribution function as 



£np(0) 



exp(l/ NP P 2 (cos6>)) 



J_j d(cos0)exp([/Np.P2(cos0)) 



(11) 



where ?7np is a variational parameter for the nanopar- 
ticles. The orientational order parameter S^p of the 
nanoparticles can be defined by analogy with the liquid- 
crystal order parameter as 

f\ d(cos 6)P 2 (cos 6) exp{U NP P 2 (cos (9)) 

Snp = -j ■ (12) 

J j g?(cos 0) exp([/NP-P2(cos#)) 

Just as in the liquid-crystal case, note that J7np ranges 
from to oo, while Snp ranges from to 1. 

As we discussed in our previous paper [17] . the ferro- 
electric nanoparticles create static electric fields, which 
interact with the dielectric anisotropy of the liquid crys- 
tal. By averaging the interaction energy over the distri- 
bution functions qlc an d f?NP, we obtain 



interaction 



(13) 



In this expression, JVnp is the number of nanoparticles 
in the system, and K^p is an energetic parameter repre- 
senting the strength of the interaction. For an unscreened 
electrostatic interaction, we derived 



K 



NP 



e Aep 2 



4:ne AsP 2 R 3 



1807r(£ e) 2 -R 3 405(e e) 5 



(14) 



where p, P, and R are dipole moment, polarization, and 
radius of a nanoparticle, and e and Ae are the dielec- 
tric constant and dielectric anisotropy of the bulk liq- 
uid crystal. If the interaction is screened by counterions, 
then -Knp is somewhat reduced, but it is still substan- 
tial as long as the Debye screening length is greater than 
the nanoparticle radius. Hence, orientational order of the 
liquid-crystal molecules tends to favor orientational order 
of the nanoparticles, and vice versa. 

Whenever there is an aligning effect, there must be 
an entropic cost. By analogy with the entropic term for 
liquid-crystal molecules, the entropic penalty for aligning 
the nanoparticles is 



F NP . 

entropic 



fc B TA r NP (ln £ NP ) 
fceTiVNP UnpSnp 



(15) 



In 



J _ 1 d(cos 0) exp([/NP-P2(cos 9)) 



The total free energy for liquid-crystal molecules and 
nanoparticles is now the combination of Eqs. ©, (|13p . 
and (fT5l). 



F 



J 



N LC k B T 



. Q2 

• D LC 



NP 



3fc B T k B T 
+£4c<Slc + j^npSnp (16) 

-In [/^(cos^exp^Lc^cosfl)) 
— ^ln d(cos9) exp(L/NP-P 2 (cos6>)) . 



Note that we have normalized this free energy by the 
number of liquid-crystal molecules, not by the number 
of nanoparticles. For that reason, all the nanoparticle 
terms in Eq. (|16|) contain a factor of v = N^p/N^c, the 
ratio of the number of nanoparticles to the total number 
of liquid-crystal molecules. 

To summarize, we have derived the free energy for the 
system of ferroelectric nanoparticles suspended in a liq- 
uid crystal. The first term represents the aligning energy 
favoring orientational order of the liquid crystal, while 
the second term describes the mutual aligning interaction 
between nanoparticle order and liquid-crystal order. The 
last terms are entropic terms that give the free-energy 
penalty for any liquid-crystal or nanoparticle order. The 
free energy is a function of two variational parameters, 
C/lc an d C^np, and we formulate our problem as mini- 
mization over those quantities. Once we find them, we 
can calculate the order parameters Slc an d ^np using 
Eqs. © and 



III. TRANSITION TEMPERATURE 

Experiments show a substantial increase in the 
isotropic-nematic transition temperature for liquid crys- 
tals doped with ferroelectric nanoparticles. In order to 
understand this phenomenon and predict how to enhance 
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it further, we investigate the isotropic-nematic transition 
using the free energy of Eq. (flU)) . 

Two distinct limiting cases of this transition are pos- 
sible. If the nanoparticle order is small, then all of the 
integrals in Eq. (|16p can be expanded in Taylor series for 
small Ulc and U^p- The expressions for the order pa- 
rameters Slc and Snp from Eqs. © and Eq. (fT2"j) can 
also be expanded in power series in J7lc and C/np- Hence, 
the free energy can be expressed as a series in ?7lc and 
JTnp, or equivalently as a series in Slc and Snp- After 
some algebraic transformations, we obtain 



N LC k B T 



= const 



J 



3fc B T 



-uS, 



NP 



knT 



SnpSlc 



(17) 



This expression is exactly the Landau free energy as a 
series in the order parameters, as discussed in our earlier 
paper 17]. To find the isotropic-nematic transition, we 
first minimize over Snp to obtain 



c Knp c 



(18) 



We then substitute this value into the free energy series 
to obtain 



F 



N LC k B T 



const - 



J 



" K NP 



2 3k B T 10(fc B T) 2 



Stc + - 



(19) 

The change in the coefficient of S 2 C shows that the 
isotropic-nematic transition temperature is shifted up- 
ward by 



AT, 



NI 



" K NP 

25fc2r NI 



(20a) 



In the notation of the previous paper, this shift can be 
written as 



NI 



ncj) N pR 3 ( 2AeP 2 ' 1 



(20b) 



where plc is the number of liquid-crystal molecules per 
unit volume and </>np = |7ri? 3 pLC jy is the volume fraction 
of nanoparticles. 

Note that the power-series approximation works well 
as long as the energetic parameter K^p is small com- 
pared with bk B T. In that case the nanoparticle order 
parameter Snp is small compared with Slc , which is ap- 
proximately 0.429 just below the isotropic-nematic tran- 
sition. However, the approximation breaks down if K^p 
becomes large compared with 5fc B T, so that Snp is large 
compared with Slc- In the latter case, the prediction for 
Snp would be greater than 0.429 on the nematic side of 
the transition. It might even be greater than 1, which 
would be unphysical. This unphysical prediction arises 
because the power-series expansion cannot take account 



of the saturation of the order parameters at low tempera- 
tures. Hence, for large -Knp we must consider a different 
limiting case. 

In the limit of large i^NP, the nanoparticle order is 
large; i.e. the variational parameter £/np approaches in- 
finity and the order parameter Snp approaches 1. In that 
case we can approximate Eq. (I12p to obtain 



Snp = 1 — 



1 



U- 



XP 



(21) 



We can then put this approximation into the free energy 
of Eq. (|16p , expand the nanoparticle entropic integral for 
large £/np, and minimize the resulting free energy over 
Unp- This calculation gives 



JJnp 
Snp 



-Knp 
k B T 



1 



Slc, 
k B T 
-^npSlc 



(22a) 
(22b) 



Note that this calculation is self-consistent, showing large 
nanoparticle order when A'np ~> k B T. Using Eqs. (|22|) . 
we obtain the approximate free energy of the nematic 
phase 



J 



N LC k B T 



3k B T 



Stc + ^lcSlc 



(23) 



-In 

^vKkp 
k B T 



/_! d(cos 9) exp(C/LC-P2(cos 9)) 
( 3XnpSlc \ 



Slc 



' In 



V 2fc B T 7 



This free energy is equivalent to the classical Maier-Saupc 
free energy of Eq. (J6]), except for the last two terms, 
which represent the energy and entropy of well-ordered 
nanoparticles interacting with the liquid crystal. These 
terms are proportional to the nanoparticle concentration 
v = Nnp/Nlc, which is small. These terms shift the ne- 
matic free energy, and hence shift the isotropic-nematic 
transition temperature. To find the value of the shift, we 
must minimize the free energy. 

To minimize the free energy, we use perturbation the- 
ory. For this calculation, we define the parameters 



Ulc = Ulr 



AUt. 



'LC T ^LC, 

Tm = T NI + ATni, 



(24a) 
(24b) 



where U^ c and T^j are the known results from the clas- 
sical Maier-Saupe free energy, given in Eqs. (J?} and ©, 
and AJ/lc and ATni are perturbations due to addition 
of ferroelectric nanoparticles. For low nanoparticle con- 
centrations, these perturbations should both be of order 
v. We now expand the free energy to lowest order in 
these pertubations, minimize over A£/lCj and solve for 
ATni such that the isotropic and nematic free energies 
are equal. The resulting shift in the transition tempera- 
ture is 



AT, 



NI 



1.03- 



vK 



NP 



k B 



1.03 



'NP 



AeP 2 



135k B pLC£o£ 2 



(25) 
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Comparing Eqs. (|20|) and (|25|) . we can see that there 
are two regimes for the shift in the transition tempera- 
ture. For small interaction K^p (i.e. the Landau regime), 
the shift ATni increases as K^ P , but for large K^p, it 
increases more slowly as K^p. In both cases it is pro- 
portional to the nanoparticle concentration v. Equiva- 
lently, if we work at fixed nanoparticle volume fraction 
0np , our theory predicts that ATni will increase with the 
nanoparticle material polarization P and radius R 3 in 
the weak-interaction regime, but it will only increase as 
P 2 and will be independent of R in the strong-interaction 
regime. (It will be independent of R as long as the par- 
ticles are small enough so that they do not distort the 
liquid-crystal alignment.) 

Our predictions for ATni can be compared with the 
previous predictions of Li et al. 13]. They calculated 
that ATni should increase as the volume fraction 0np 
and as the polarization P 2 , and should be independent 
of the radius R. These predictions for the scaling agree 
with our predictions for the strong-interaction regime (al- 
though not for the weak-interaction regime). We believe 
that this agreement is just a coincidence, because the the- 
ories are quite different. One way to see the difference 
is through the dependence on dielectric anisotropy Ae: 
they predict that ATni should scale as (Ae) 2 , but we cal- 
culate that it should scale linearly with Ae in the strong- 
interaction regime. This difference arises because their 
model considers one liquid-crystal molecule interacting 
through the dielectric anisotropy Ae with one nanoparti- 
cle, which then interacts through Ae with another liquid- 
crystal molecule, thus giving an effective liquid-crystal in- 
teraction proportional to (Ae) 2 . By comparison, in the 
strong-interaction regime our model considers the direct 
influence of well-ordered nanoparticles on the liquid crys- 
tal, and hence has only one power of Ae. 

For a numerical estimate, we use typical experimental 
values of the parameters 0np = 0.5%, P — 0.26 Cm~ 2 , 
R = 35 nm, p LC = 2.4 x 10 27 m" 3 , fc B = 1.38 x 10~ 23 
JK" 1 , e = 8.85 x 10" 12 C^N^uT 2 , and Ae w e « 10. 
Those parameters imply v = 1 x 10~ 8 , K^p = 1 x 10~ 15 
J, and hence K^p/{kp,T) = 2x 10 5 , so the system is def- 
initely in the strong-interaction regime. Our prediction 
for the shift in transition temperature is then 



AT 



XI 



1 K. 



(26) 



This value is consistent with the order of magnitude that 
is observed in experiments. Note that in this prediction 
we are using the bulk polarization of the ferroelectric ma- 
terial BaTi03, which is P = 0.26 Cm~ 2 . In this respect, 
our current estimate is different from our previous pa- 
per pj}, where we assumed P — 0.04 Cm -2 because of 
an understanding that the bulk polarization is reduced by 
surface effects in nanoparticles. The issue of estimating 
the polarization of nanoparticles is subtle, as discussed 
inRef. [Tjj|. 

As a final point about the phase diagram, we should 
mention that the model defined by the free energy (|16[) 
can exhibit one additional phase, between isotropic and 



nematic, which occurs if the parameter -Knp is suffi- 
ciently large. In this intermediate phase, the nanoparti- 
cles have substantial orientational order (with Snp com- 
parable to the Maier-Saupe order parameter of 0.429), 
but the liquid crystal has only very slight orientational 
order (with Slc of order vK-^p / {kftT)). For that reason, 
we might call it a "semi-nematic" phase. It is a pertur- 
bation on the pure liquid crystal's isotropic phase, not on 
the nematic phase. The semi-nematic phase is probably 
an artifact of the mean-field theory used here. It can only 
exist because the very slight order of the liquid crystal 
mediates an aligning interaction between the nanoparti- 
cles. This slight orientational order is unlikely to persist 
when one includes fluctuations in the liquid crystal. 



IV. KERR EFFECT 

Apart from the phase diagram, another important is- 
sue is the response of a liquid crystal to an applied elec- 
tric field. In the isotropic phase, an applied field E in- 
duces orientational order proportional to E 2 , known as 
the Kerr effect. In most pure liquid crystals, the Kerr 
effect is quite small, and can only be observed for very 
large fields. However, in our previous paper, we predicted 
that ferroelectric nanoparticles can enhance the Kerr ef- 
fect by several orders of magnitude. We would like to as- 
sess how this prediction is modified by the Maier-Saupe 
theory presented here. 

In the presence of an electric field, ferroelectric 
nanoparticles will have polar order along the field; i.e. the 
orientational distribution function will no longer have a 
symmetry between the directions +11 and — n. Hence, we 
must change the nanoparticle distribution of Eq. (|11[) to 



0np(0) = 71 



!7f p Pi(cos e) + (7j ,p P 2 (cos 6) 



J_ 1 d(cos9)e u ? rp ^ cose ') +U 2 rp ^ cose ') 



(27) 



Here, Uf p and are two variational parameters, 

which act as effective fields on the polar and nematic or- 
der of the nanoparticle distribution function, as described 
by the Legendre polynomials Pi (cos 9)) and T2(cos#)), 
respectively. They generate polar and nematic order pa- 
rameters, defined as 



NP 



S 



NP 



f\ d(cos 8)P X (cos 9)e u ^ p ^ ( cos e)+urP^os e) 
J* d(cos 9)e u i Ppl ( cos e )+^2 NP ^ (cos J) ' 

(28a) 

j\d{cos9)P 2 {cos9)e u ^ p ^ cos ^+ u ^ p ^ cos9 '> 



J_ 1 d(cos9)e u i Pp ^ cose '>+ u 2 Fp ^ f >) 



(28b) 



We still assume that the liquid-crystal distribution func- 
tion is purely nematic, not polar, as given by Eq. ([2]l. 
The applied electric field E adds two contributions to 
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Saturated LC order 




io- e I 

10" 5 1CT 4 10 s tO"* 10" 1 10° 10 1 10' 10 3 10* 10* 

Electric Field (V/|im) 

FIG. 2: (Color online) Four regimes of the Kerr effect, de- 
rived from a numerical minimization of Eq. (I30p with typical 
experimental parameters. A log-log scale is used to show all 
the regimes on a single plot. 



and (|20| . In this expression, the first term is the conven- 
tional Kerr effect without nanoparticles, and the second 
term is an additional contribution due to the aligning ef- 
fect of the nanoparticles. Note that both terms are pro- 
portional to E 2 . With the numerical estimates presented 
above, the second term is several orders of magnitude 
larger than the first, and hence the nanoparticles greatly 
increase the Kerr effect in this regime. 

(b) For larger field, in the regime k B T/p < E < 
{vK^ppLc/^eoAe)] 1 / 2 , the nanoparticle order parame- 
ters -Mnp and Snp saturate near the maximum value of 
1. In that case, we can no longer expand the free energy 
as a power series in the nanoparticle parameter, but we 
can still expand it in the liquid-crystal parameter. Min- 
imizing the free energy then gives 

^= i5 M r-r.) (^ + ^4 (32) 



the energy of the system 
e Q Ae 



afield 



E 2 Slc-Nlc 



pEM^pNjqp. 



(29) 



Here, the first term is the interaction of the field with the 
dielectric anisotropy of the liquid crystal, and the second 
term is the interaction with the dipole moments of the 
nanoparticles. With these energetic terms, together with 
the entropy of the distribution function, the free energy 
becomes 



F 



N LC k B T 

e AeE 2 



J Q 2 _ vK ^P „ 



lc-Snp 



S 



3k B Tpi,c ~ LC k B T 
+U LC S LC + vU? F M NP + i/[/ 2 NP S NP 



vpE 



NP 



(30) 



In 



-is In 



f\d(cos6)e u ^ p ^ cos ^ 
J_, d(cos9)( 



The next step is to minimize this free energy over all 
three variational parameters Ulc, Uf p , and . For 
this minimization, there are four distinct regimes of elec- 
tric field, as indicated in Fig. 2. 

(a) If the field is sufficiently small, E < k B T/p, it in- 
duces only slight order in the liquid-crystal and nanopar- 
ticle distributions. In that case, we can expand the free 
energy as a power series in all the variational parameters. 
This expansion is exactly the Landau theory presented in 
our previous paper [l7j . We can then minimize the free 
energy over all the variational parameters to obtain 



Slc = 



E 2 



15fc B (T-T* ) V PLC 5(*bT) 



e Ae vK^pp 2 



(31) 



where T£ d = T* + AT NI is the limit of supercooling of 
the nanoparticle-doped liquid crystal, combining Eqs. @ 



Once again, the first term is the conventional Kerr ef- 
fect without nanoparticles, and the second term is the 
additional contribution from the nanoparticles, but now 
the second term is independent of electric field. The sec- 
ond term is still much larger than the first, and hence 
the Kerr effect is approximately constant with respect to 
field in this regime. 

(c) For even larger field, [^ N ppLc/(eoAe)] 1/2 < E < 
[k B (T - T*)p LC /(e Ae)] 1 / 2 , the order parameter Slc is 
still given by Eq. (|3"2"1) . but now the first term becomes 
larger than the second. In this regime, Slc again in- 
creases as E 2 . It is similar to the conventional liquid- 
crystal Kerr effect, but with an extra constant contribu- 
tion from the nanoparticles. 

(d) For the largest field, [k B (T - T*)p hc / {eqAe)] 1 / 2 < 
E, the order parameter Slc saturates at the maximum 
value of 1. 

To get a full picture of the behavior through all these 
regimes, we minimize the free energy of Eq. (|30[) numer- 
ically, using the typical experimental parameters listed 
at the end of Sec. III. The results of this calculation are 
shown by the black line in Fig. 2. By comparison, the red 
line shows the limiting case of regime (a), and the green 
line shows the approximation for regimes (b-d). We see 
that the numerical solution overlaps the limiting cases 
and connects them. 

Note that the low-field regime (a) is the regime where 
Landau theory is valid, and it is where the nanoparticles 
give the greatest enhancement of the conventional Kerr 
effect. However, this regime will be difficult to observe in 
experiments, because the induced order parameter Slc 
is so small, on the order of 10~ 4 . Typical optical ex- 
periments can only detect a birefringence corresponding 
to Slc 011 the order of 10~ 2 , which does not occur until 
regime (c), which is closer to the conventional Kerr effect. 
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V. CONCLUSIONS 

In our previous paper [I?} , we developed a Landau the- 
ory for the statistical mechanics of ferroelectric nanopar- 
ticles suspended in liquid crystals. This theory differs 
from other models by considering the orientational dis- 
tribution function of the nanoparticles as well the liq- 
uid crystal. It shows a coupling between the nanoparti- 
cle order and the liquid crystal order, which leads to an 
increase in the isotropic-nematic transition temperature 
and in the Kerr effect. In the current paper, we consider 
the same physical concept, but we improve the mathe- 
matical treatment by using a Maier-Saupe-type theory. 
This theory reduces to the previous Landau theory in the 
limit of weak interactions (for the isotropic-nematic tran- 
sition) or weak electric helds (for the Kerr effect). How- 
ever, it changes the results in the opposite limit, when 
the order parameters begin to saturate. For that reason, 
the new theory should make more accurate predictions 



for experiments. 

In general, the concept of coupled orientational dis- 
tribution functions should be useful for many other sys- 
tems beside ferroelectric nanoparticles in liquid crystals. 
For example, it applies to any type of nonspherical col- 
loidal particles, such as carbon nanotubes, in a liquid- 
crystal solvent. It also applies to two distinct species of 
nonspherical colloids suspended in an isotropic solvent, 
which could have a coupled ordering transition. Such 
systems would provide further opportunities to investi- 
gate the theory presented here. 
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